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We test application of the maximum entropy method to decompose the states contributing to the unstable g 
correlation function through the spectral function in the four dimensional 0(4) (/>^ theory. Reliable results are 
obtained for the a mass and two-particle tttt state energy using only the a correlation function. We also find that 
the property of the a particle is different between the unstable {ma/rriTr > 2) and stable (mo-Zm^ < 2) cases. 
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1. Introduction 

The decay of particles, e.g., p meson decay is 
not understood well on the lattice. One of the 
difficulties is that an unstable particle and its in- 
termediate states, such as the p meson and the 
TTTT state in the isospin-1 channel, have the same 
quantum numbers, leading to a multi-exponential 
form of the unstable particle correlation function. 
To obtain the energies of the states we have to de- 
compose the states in the correlation function. 

In a finite volume it is possible to decom- 
pose the states with the spectral function, which 
can be numerically reconstructed from the cor- 
relation function with the maximum entropy 
method(MEM). MEM has been recently applied 
to meson correlation functions in quenched lattice 
QCD [1]. Here we explore application of MEM to 
calculate the energies of states in unstable parti- 
cle systems. 

2. Model and parameters 

Since the calculation of the unstable particle 
is difficult in QCD, we use the four dimensional 
0(4) (f>* theory in this work. In this theory the 
a particle and two-pion state in the isospin-0 
channel have the same quantum numbers. These 
states overlap in the a correlation function. 

We calculate the correlation functions for a 
and TTTT with zero momentum for the two cases, 
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one where a is unstable {ma/mT^ w 3.7) and the 
other stable {ma/rriT^ « 1.8). Several spatial lat- 
tice sizes in the range 10^ — 28^ are employed 
to study the volume dependences of the spec- 
tral functions, while the temporal lattice size is 
fixed to 64. We perform 0.6 to 1.2 million itera- 
tions per simulation point and lattice size. The 
method of simulation and measurements of the 
correlation functions follow ref. [2] . We calculate 
((0(t) - 0(r + l))O(0)) to subtract the vacuum 
contribution where O denotes a or tttt operator. 

3. Results 

3.1. Spectral function data 

In Fig.l we illustrate the a and tttt correlation 
functions at two volumes in the unstable case. 
For smaller volume the slopes of the two correla- 
tion functions agree, while for larger volume the 
slopes are different and the a correlation function 
is of multi-exponential form. This means that the 
overlap of the tttt state dominates the a correla- 
tion function for small volume, while it decreases 
as volume increases. 

The spectral function reconstructed from the 
correlation function at various volumes in the un- 
stable case is presented in Fig. 2. The parame- 
ters for MEM are compiled in ref. [3]. We find 
the spectral functions to be a sum of (5-function 
peaks as expected. 

The first peak corresponds to the tttt state with 
zero momentum and the second peak to the a 
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Figure 1. Correlation functions for the a and tttt 
in the unstable case. 

state. As seen from the peak height, the over- 
lap of the TTTT state in the a spectral function de- 
creases with increasing volume. This explains the 
difference of behavior of the a c;orrelation function 
in Fig. 1. The a peak in the tttt spectral function 
also decreases as the volume increases. 

In the a spectral function the downward arrow 
marks the position where the tttt state with rel- 
ative momentum p = ^-k/L should appear. The 
state is not observed, showing that the overlap is 
too small in these small volumes. 

3.2. Sigma mass and two-pion energy 

In the unstable case, the a mass and tttt en- 
ergy are obtained from the positions of the peaks 
in both the a and tttt spectral functions. We 
compare results from this analysis with those ob- 
tained with the diagonalization of the correlation 
function matrix [4]. The results are shown in 
Fig. 3; MEM(c7) and MEM(7r7r) denote data ob- 
tained from the u and tttt spectral functions, re- 
spectively. For larger volumes the a mass from 
MEM(7r7r) and the tttt energy from MEM(c7) suf- 
fer from large errors or do not appear, because 
the overlaps of the states are small for large vol- 
umes. We find the a mass from MEM to be rea- 
sonably consistent with that from the diagonal- 
ization. For the tttt energy all results agree well. 

In the stable case we also obtain the a mass 
and TTTT energy with MEM and diagonalization. 
The results are presented in Fig. 4. The two set 
of data are reasonably consistent. 

Comparing Fig. 3 with Fig. 4, we find that 
there is an essential difference in the volume de- 
pendence of the a mass and tttt energy between 
the unstable and stable cases. In the unstable 
case the a mass decreases as the volume increases 
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Figure 2. Spectral functions for the a and tttt in 
the unstable case. 

while it increases in the stable case. The tttt en- 
ergy also shows a different volume dependence for 
the two cases. 

3.3. Volume dependence of 

In order to understand the volume depen- 
dence of cr mass, we fit our data from di- 
agonalization assuming the perturbative form 
in a finite volume given by m„{L) = ma + 
gii{Amcr{L) — Amcr(oo)) where Amcr{L) = 



1 

with 

£1 = 1- 
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+ 2C^{p)+6Ca{p) 



Wc(p)(niJ-VK^(p)) ' 

and Wa{p) = 2yjm'i+f-, pi = 2sin(pi/2), pi = 
2-Kni/L. Here nia- is a mass in the infinite vol- 
ume and gu is the renormalized coupling, which 
are the fit parameters. We employ simulation re- 
sults for rriTT obtained from the pion correlation 
function at L = 28. The fit results are shown by 
dashed lines in Fig. 3 and Fig. 4; they agree quite 
well with the simulation results. We note that 
the difference in volume dependence is caused by 
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Figure 3. Energies in the unstable case. 

the sign of C^{p) between the unstable and stable 
cases. 

3.4. Volume dependence of Et^tt 

To understand the volume dependence of tttt 
energy, we use the scattering length oq that 
relates to the difference of Et^tt from 2m^ 
through Liischer's formula [5] given by Et^tt — 

= + +C2#) , with Ci = 

-2.837297, = 6.375183. We employ the tttt en- 
ergy obtained with diagonalization. The results 
are presented in Table 1. The scattering length 
exhibits an opposite sign between the two cases. 

We calculate the perturbative scattering length 
from the perturbative phase shift [2] and find 

"0 = 9fom4 fl^'- 4 ' where R = ma/m^. The 
renormalized coupling gjt is already determined 
in the a mass fit. Alternatively one may use 
the perturbative relation gn = 3Z„{m'^ — m.'^)/v^ 
where m,r, and v are taken from the simulation 
results for _L = 28. We employ the a mass deter- 
mined from the a mass fit for nia-- We compile ao 
and gn from simulations and from perturbation 
theory in Table 1. 

We find results from the simulation to be rea- 
sonably consistent with the perturbative results. 
The sign of the perturbative scattering length de- 
pends only on whether the a particle is unstable 
or stable. Prom these signs the tttt energy in a 
finite volume is expected to increase or decrease 
in the unstable or stable cases, respectively. The 
behavior of the tttt energy also agrees with our 
simulation results. 



0.8 

0.7 

0.6 

0.5 

0.4 

0.3 
0.5 

0.4 

0.3 

0.2 



X diagonalization 
• MEM(c) 
O MEM()iJt) 



X3 




Figure 4. Energies in the stable case. 
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9r 



ao 



Qr 



(S) diag. 0.289(9) -2.49(19) 

(P) 5fl(fit) 0.383(3) 14(1) -1.98(23) 9(1) 

(P) fffl(def.) 0.496(1) 18.9(2) -3.62(11) 17.8(1) 

Table 1 

Scattering lengths obtained from simulation (S) 
and perturbative formula (P). 

4. Conclusion 

We have demonstrated that the maximum en- 
tropy method is a simple technique for decompos- 
ing the states in the unstable particle correlation 
function and determining their energies. In future 
it will be interesting and challenging to apply this 
method to studies of the p meson decay and the 
a particle in QCD. 
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